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Abstract We present some exact results on the behavior of Branching and 
Annihilating Random Walks, both in the Directed Percolation and Parity 
Conserving universality classes. Contrary to usual perturbation theory, we 
perform an expansion in the branching rate around the non trivial Pure Anni- 
hilation model, whose correlation and response function we compute exactly. 
With this, the non-universal threshold value for having a phase transition in 
the simplest system belonging to the Directed Percolation universality class 
is found to coincide with previous Non Perturbative Renormalization Group 
approximate results. We also show that the Parity Conserving universality 
class has an unexpected RG fixed point structure, with a PA fixed point 
which is unstable in all dimensions of physical interest. 

Keywords Reaction-diffusion • Exact results • Phase transitions • Non- 
equilibrium • Field theory 



1 Introduction 

The study of critical behavior in out of equilibrium systems has been a very 
active topic in statistical mechanics during the last decades [TJ[2J|3]. As in 
equilibrium, fluctuations and correlations become large in systems close to 
a continuous phase transition, leading to divergences in quantities such as 
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correlation time and length, and to emergent phenomena classifiable (as in 
equilibrium) in different universality classes. 

Renormalization Group (RG) methods have been employed since its de- 
velopment to the study of critical dynamics pQ . These are well suited for the 
task, given that the diverging correlation length in second order phase transi- 
tions signals the emergence of scale free behavior, whereas the RG approach 
focus on how systems change under scale transformations. 

When the system in question is far from equilibrium, as for example 
when it is not constrained by detailed balance, it usually exhibits a much 
richer variety of phenomena than their counterparts close to equilibrium. In 
particular, usual RG techniques have to be adapted, given that one does 
not know the explicit probability distribution, in the way one knows the 
Boltzmann-Gibbs distribution when in equilibrium. 

Perturbative RG has been used in the study of second order out of equi- 
librium phase transitions, although not with the same level of success as in 
the study of equilibrium phase transitions. This is due in part to the absence 
of high-order perturbative results, as opposed to equilibrium problems, but 
also to deeper physical issues. For example, for most of these systems, upper 
critical dimensions happen to be generally far from the dimensions of physical 
interest, difficulting the usual e-expansion. Moreover, these models generally 
lack a lower critical dimension or (generally speaking) an exactly solvable 
low d model. Finally, out of equilibrium systems also tend to be more prone 
to show genuinely non-perturbative behavior, such as large couplings [5, 5 > 6 . 

In this work we study some properties of phase transitions occurring in 
what arc known as Branching and Annihilating Random Walks (BARW) [3J 
I7U81I91IT0] . that is, systems composed of particles of a single species A, which 
diffuse in a d-dimensional space, and that can suffer both annihilation and 
branching (i.e. offspring creation) processes, with different rates. From these 
competing processes typically emerges, at long times, a stationary state which 
can either be in an active or an absorbing phase, with the absorbing phase 
corresponding to a no-particles, no-fluctuations state. This implies in partic- 
ular the absence of detailed balance in these systems. The transition between 
both phases, which can take place depending on the microscopic rates, is typ- 
ically of a continuous type. BARW are not only of direct physical interest, 
but also present a relatively simple class of out of equilibrium systems, and 
have shown to be very useful for the study of the role of fluctuations in out 
of equilibrium statistical physics [3]. 

Due to universality, it is in general enough to consider the simplest pos- 
sible reactions, as for example 

2A 4 0, A ^ 2A (1) 

or 

2A \ 0, A A 3A (2) 

(when A —> 2A does not exist in the second case) as these reactions are 
the most relevant in the RG sense. At the mean field level (that is, for the 
classical rate equations) no phase transition is found. That is, fluctuations 
are here responsible for the very existence of a phase transition. This is in 
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stark contrast with most other known phase transitions, where mean field 
results predict in general the presence or absence of a phase transition, even 
when they are unable to yield accurate phase diagrams or critical exponents. 
In the case of BARW, the mean field result shows the need for taking into 
account statistical fluctuations, a task for which one expects RG methods to 
be particularly effective. Notice that in the definition of BARW we exclude 
explicitly the reaction A — > 0. If such reaction is present, the phase diagram 
is qualitatively well described at a mean field level. Notice also that in this 
work we will not be dealing with the Pair Contact Processes with Diffusion 
(PCPD) universality class [TTlll2j . which can be seen, in terms of BARW, as 
systems whose reactions involve always at least two particles. 

BARW can be classified into sub-classes [3M], and in this work we will 
concentrate on the simplest two, which depend on the presence or absence of 
a symmetry conserving the parity of the number of particles. If no such sym- 
metry exists, as is the case of (UJ, it has been shown that the BARW system 
belongs to the Directed Percolation (DP) universality class [Tj5] whenever 
a second order phase transition takes place. When only reactions preserving 
the parity of the number of particles (e.g. ©) are present, an additional sym- 
metry appears, changing the universal properties of the system. If a phase 
transition takes place in this case, it is known to be in the Parity Conserv- 
ing (PC) universality class [3] (also more properly known as Generalized 
Voter universality class |12j). From now on we will refer to these systems as 
BARW-DP and BARW-PC respectively. 

There exist various known results about BARW, and even some exact 
results for low-order vertices [Tl] or for special BARW systems which do not 
present phase transitions |15j . Within perturbative RG, a phase transition 
for the simplest BARW-DP system (the one consisting of the reactions ([1])) 
is found for space dimensions d < 2 only [8 . This improves the mean field 
result, but still contradicts Monte-Carlo and Non Perturbative Renormaliza- 
tion Group (NPRG) results, which observe a phase transition for any d [TBI 
117] . This difficulty of the perturbative approach may have to do with the 
fact that for d > 2 the transition occurs for values of the annihilation rate A 
which are large, and thus out of reach of a perturbative analysis performed 
around the reaction-less Gaussian fixed point. 

As for BARW-PC, the perturbative studies of [5] showed the existence of a 
new universality class different from DP, and of a new fixed point for d smaller 
than a new critical dimension d c ~ 4/3. This behavior had already been 
predicted as a consequence of the additional symmetry [7J, although some 
early studies [18| were not conclusive with respect to this new universality 
class. In the NPRG context there have also been studies of BARW-PC [19] . 
which seem to confirm the existence of a new fixed point for d < d c ~ 4/3. 
In both scenarios the appearance of this new fixed point is associated with 
a change of stability of the Pure Annihilation (PA) fixed point: for d > 
d c , a is relevant and only an active phase exists, whereas for d < d c , a is 
irrelevant and an absorbing phase with the properties of PA at long distances 
is expected for small a . 

In this work (i) we obtain exact and closed equations for all response 
functions of the PA model and (ii) we show how to perform an expansion 
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in a around this model. Since our approach is valid for any value of the 
annihilation rate A, we obtain exact results at small a. 

For BARW-DP we show that an active to absorbing phase transition 
exists in all dimensions d and we compute the non-universal threshold values 
Xth{d) above which it occurs for two specific microscopic realizations of the 
system. 

For BARW-PC we show, in disagreement with both the 2-loop perturba- 
tive and Local Potential Approximation (LPA) results, that the stability of 
the PA fixed point does not change between one and two dimensions. This 
contradicts the existing scenarios explaining the existence of a phase transi- 
tion in d = 1. We propose an alternative scenario that reconciles all existing 
results. It is based in the appearance in a dimension d €]1,2] of two fixed 
points that split apart as d is decreased. The fixed point with the small- 
est fixed point value for the branching rate should be fully attractive, while 
the other one should have a repulsive direction associated with the active to 
absorbing phase transition belonging to the PC universality class. 

Some of these results have been recently presented in an abridged form 
in [20]. 

The paper is organized as follows. In Section [5] we give a quick overview 
of the application of field theory to the study of out of equilibrium statis- 
tical systems. The interested reader is nonetheless strongly advised to read 
more general reviews of the these methods [3101121] . In Section [3] we show a 
method to find all generalized response functions in the steady state of the 
simple reaction diffusion system corresponding to PA. In Sections [4] and [5] we 
propose an expansion around this solution, in order to analyze BARWs in 
both universality classes, and to answer some specific questions concerning 
their phase diagrams. We have decided, in order to make the proofs simple, 
to use along the main part of the article derivations based on resummations 
of perturbative series. For completeness, however, we give non perturbative 
proofs (beyond an all-loop order analysis) of our results in the appendices, 
as well as presenting some other technical details. 



2 Field theory for BARW 

There are many known methods in the literature for the mapping of out of 
equilibrium problems onto field theories [2TII221I23"] . In the case of reaction- 
diffusion processes, a field theory can be constructed in a standard way by 
using the Doi-Peliti formalism [23], the idea of which is to re-express the 
Master Equation for the occupation probabilities in a lattice system using 
creation and annihilation operators in an abstract Fock space, followed by 
a coherent-state path integral representation, and (optionally) the use of a 
continuum limit for the lattice. As a result of this procedure, one obtains a 
functional integral (the so-called generating functional) 

Z[J, J] = J V(f>V(j) exp (-S[(/>, </>] + J J<f>+ J^j (3) 
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with an appropriate action S[4>,(j)], which captures exactly the microscopic 
reactions. Here we have introduced the notation, to be used throughout 

x = (x, t) and p — (p, v) (4) 

where the last convention will be used in Fourier space. We also introduce 
some notation for the integrals 

The time-dependent statistical correlation and response functions can 
then be computed from iJ, by functional derivation w.r.t. to the sources J 
and J. In this context, the expected value of the field <fr(x) is associated with 
the local density of A particles, and the response field (f)(x) allows for the 
computation of response functions. 

For general processes of the type A (m 4- \)A and kA — ^» 0, with 
diffusion constant D, this procedure yields (ignoring initial conditions, which 
play no role in the long time stationary state, see for example [9]) 

S[M]= [ U(x)(d t -DV 2 )4>(x)-\ k (l-$(x) k )4(x) k 

+ a m (l-4>(x) m )4>(x) ( f>(x)) (6) 

Diffusion is responsible for the kinetic part (corresponding to Brownian 
motion). Reactions give rise to interaction terms in the potential-like part of 
the action. A perturbative expansion can be set around the exactly solvable 
Gaussian part of the action in the usual way |24j . One can perform a pertur- 
bative expansion [5] to approximate the correlation and response functions 
of the theory at any desired order in a m and e = d c — d, with d c the 
upper critical dimension, above which mean field results are expected to give 
a good description of the universal properties of the system. As stated in 
the introduction, for the purposes of the present work we will concentrate 
on reactions involving a minimal number of particles. These are enough to 
characterize the universal properties of these systems, and, for non-universal 
properties, they can be seen as the simplest examples. 

The connected correlation and response functions of a theory will be 
written 

G^ n ' m \x 1 ,...,x n ,xi, ...,%) = {(t>(xx) . . . <t>(x n )<f>(xi) . . . <f>(x m )) c (7) 

which are generated by taking derivatives of the logarithm of the generating 
functional log Z [J, J] . These can be obtained in a perturbative series using 
connected Feynman diagrams. In this work we will mostly work with the 
vertex functions r^ n ' m \ the amputated 1PI functions of the theory, which 
include all the information coming from fluctuations in the system. The gen- 
erating functional -T[(</>), (</>)] for the _T(™' m ) vertices is given by the Legendre 
transform of the connected generating functional log Z [J, J]. 
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It is often convenient to perform a shift in the fields, of the form [8] 



which allows for some simplifications in the functional form of the interaction 
potential, and is needed to make the nexus between BARW-DP and Directed 
Percolation. This shift is not convenient in the BARW-PC case however, 
where it obscures the presence of the related parity conserving symmetry. 

In the case of out of equilibrium models, special care must be taken with 
respect to the causal structure of the theory. In this regard, actions such as 
(H2), given by the Doi-Peliti formalism, implicitly require the use of the Ito 
prescription, in which all quantities are evaluated with the convention that 
the Heaviside function 0(t) is zero for t = [251126) . In perturbation theory it 
is relatively easy to implement the so-called Ito prescription, as it amounts to 
force closed propagator loops to be zero [51I2T]. Non-perturbative equivalent 
results are given in [3] and in Appendix |DJ 

Diagrammatically, each (jxp propagator can be represented as a line with 
an arrow going from to (f>, and each such propagator carries a Heaviside 
function of time, expressing causality. We use in the following the diagram- 
matic convention of drawing only <p-(p propagators (that is, the function 

G(p) = [.T^ 1 ' 1 ) (— p)l ) but we include, if allowed in a given model, insertions 
of rt 2 '^ or r< ' 2 ) as vertices. 



3 Pure Annihilation 

In this Section we study the simplest case of a reaction-diffusion system, PA, 
in which the only reaction in the system is annihilation by pairs of diffusing 
particles A + with a probability rate A. Later we will use the exact 

solution for this particular system as the starting point of a perturbative 
expansion, in order to study more general BARW at small branching rates. 
It is easy to prove [9"I[T1] that this system belongs to the same universality 
class as pure coagulation, in which the only reaction is A + A — > A. In the 
following we will use the PA model but the pure coagulation case can be 
analyzed in a similar way. 

After implementing the Doi-Peliti procedure and performing a shift in 
the response fields, Eq. (O, the bare action S PA can be written [HIE] 



As said before, we only analyze the steady state where all correlation func- 
tions are zero, since the system always approaches the empty state in the long 
time limit. However, even in this state, the response functions are non triv- 
ial, and are governed in the IR (that is, for momenta and frequencies smaller 
than the scale set by A) by a non-trivial fixed point of the RG equations, for 
d < 2. In the following we speak of 'correlation functions' in a generalized 
sense, including response functions. 



<f>{x) = 1 + 4>{x) 



(8) 



S PA [4>, <f\= ($(d t - DV 2 )0 + \^{4> + 2)0 2 ) 



(9) 
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As it stands, this theory shows a certain resemblance with the standard 
<j) A scalar field theory. However, symmetry and causality properties allow for 
a greatly simplified analysis. We first show that for the PA model all r^ n,m ^ 
functions can be obtained from the r( n,n \ vertices with the same number 
of incoming and outgoing legs. This is quite clear perturbatively, but we 
give in the following a non-perturbative proof based on a Ward identity for a 
U(l) transformation. Secondly, we deduce a general identity yielding a closed 
equation for any J^™*™). It is easy to verify that the T^ 2 ' 1 ) and 

vertices thus obtained coincide with the results of (81,14,27,28!. We show in 
Appendix lAl how to compute _r( 3 ' 3 ) from our method. 



3.1 U{1) Ward identity 

Let us start by studying a generalization of PA with action S PA , where 
couplings for the cubic and quartic terms are independent. 

S PA $A] = J (4>(d t -DV 2 )<t> + A 3 ^ 2 + A 4 (^) 2 ) . (10) 

Let us consider the Ward identity (51] associated with the infinitesimal field 
transformation 

4>{x) -> (l + e)4>{x) 

4>{x) -> (1 - e)4>{x) (11) 

When A3 = this is a symmetry of the action, but the cubic term breaks 
it explicitly. We can nevertheless obtain a Ward identity associated with this 
transformation by performing (1111) as a change of variables in the expression 
for Z [J, J] , given in Eq. ([3]) : 



= e (Jt-W + Wfljj (12) 

J X 

Here the mean value (■ • • )j j is computed in the presence of the sources J 
and J. The term proportional to A3 can be written as a derivative w.r.t. A3 
of the generating functional of connected correlation functions. By Legendre 
transforming Eq. (fl2"j) . one deduces the Ward identity 

+ f(& -*%)-<> (is) 



This equation can be derived w.r.t. <f> and (f> fields and evaluated at zero field, 
yielding 

(n-m)r^ = \ 3 ^—- (14) 
CA3 

where r^ n,m ^ is a function of (xx, . . . , x n , x%, . . . , x m ). 

Since, perturbatively, 7^™'™) can only involve positive powers of A3, this 
equation shows that J T(">" 1 ) - 0(X^~ m ) when n > m, and that for PA r("- m ) 
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contains exactly (n— to) third-order bare vertices. We conclude that all I^™' n ) 
vertices can be computed directly from the action with A3 = and that 

r^ n ' m \x 1 ,...,x n ,x 1 ,...,x m ) =0 ifn<TO (15) 



which simplifies the study of this system. 

Given these results one can conclude that for any correlation function, 
the perturbative expansion in A3 is, being in fact a polynomial, exact at a 
finite order. In order to calculate the connected correlation function G tyn,m ^ 
(with n > to), one can expand the functional integral at order X^~ m : 

{x\ , . . . , x n . . . ,x m ) 

_ _ / f _ \ n—m 

Ar m (#zi) ■ • ■ 4>(x n )^{ Xl ) . . . 4>{x m ) ( / U 2 ) >c 



(16) 

=./=n.A a =n 



(using the unique decomposition of G"( n ' m ) in terms of 1PI vertices [53]) 
reducing its calculation to the knowledge of correlation functions of the A3 = 
model (which only contains J^™*™) vertices). This shows that the building 
blocks of the PA model are the vertex functions with an equal number of 
incoming and outgoing legs that can be calculated at A3 = 0. 



3.2 An identity for the r {n ' m) vertices 

We now present an identity allowing us to obtain a closed equation for any 
p(n,m)_ can kg mos t conveniently written at the diagrammatic level: any 
diagram contributing to / n ("' m ) which includes at least one loop has the 
structure shown in Fig. [T] (that is: any 1PI perturbative diagram begins with 
a 4-legs bare vertex). 




n — 2 legs 



Fig. 1 Generic form of a diagram contributing to r^ n ' m > and that involves at least 
one loop in PA. Left hand side: diagrammatic representation of a generic r (n ' m ^ 
vertex. Right hand side: general structure for such vertices in PA, the black blob is a 
connected and amputated Green function that has to comply with some requisites, 
see text. 



9 



The black blob denotes a sub-diagram that is constrained by the condition 
that the full diagram must be 1PI. In particular, it means that this sub- 
diagram must be connected (and with amputated external legs). Now, any 
connected diagram with n incoming and m outgoing legs has a unique tree 
decomposition in terms of 1PI sub-diagrams having at most these numbers 
of legs. By summing all possible diagrams and permutations compatible with 
the 1PI structure of the full diagram, we obtain a closed equation that re- 
lates any r^ n,m ' with vertices I^ s '') with a lower number of legs. A non- 
perturbative proof (not based on an all-order analysis) of this general prop- 
erty is given in Appendix [Cj by using NPRG techniques (to which we give 
an introduction in Appendix FBI) . 

Notice that, as explained in the previous Section, 7^™'™) vertices can be 
calculated at A3 = 0. Now, for A3 = the U(l) transformation ITT]) is a sym- 
metry of the action, and the fields <f> and (f> play a symmetric role. Accordingly, 
this same construction can be performed singularizing two outgoing legs in 
the case of J^™'™) vertices. 

In order to be concrete, let us analyze the identity given in Fig. [T] for 
the simplest vertices. For i"^ 1 ' 1 ) this gives a well-known non-renormalization 
property: there is no correction to 2" 1 ' 1 ' 1 ) in PA. This is due to the fact that 
no diagram such as the one presented in Fig. [T] can be drawn with a single 
incoming leg. This no-field-renormalization condition implies that the critical 
exponents r\ and z have their mean field values, r\ — and z = 2. Concerning 
r^ 2,2 \ the result is less trivial. Given that there are only two incoming legs 
and that in a theory without cubic vertices all connected diagrams with four 
external legs are 1PI, one arrives at the closed equation (see Fig. [2]) 




Fig. 2 Closed equation for r (2,2) in PA. 



which reads: 

rW(p 1 ,p 2 ,p 1 ,p 2 ) = 4A 4 - 2A 4 f G{q) 

Jq 

X G(pi + P 2 - q)r {2 ^ (q, p 1 +p 2 -q,p 1 ,p 2 ) (17) 
whose solution is of the form (see Appendix [E] for a proof) 

rW(p u p 2 ,pi,p 2 ) = Ah( Pl + P2 ) (is) 

By substituting in TO 

we find 

UiP) = l + 2M fl ^)G(p-,) (19) 
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For i^ 2 ' 1 ) the identity in Fig. [T] becomes that of Fig. [3l which can be 
written as 

rW(p 1 ,p 2 ,p)=2\ 3 -\ 4 f G(q) 

x G( Pl + P2 - q)rW ( Pl + P2 ~q, q, p) (20) 

We can show (see Appendix IB"]) that this implies that I^ 2,1 ) (pi,p2,p) depends 
only on p. We thus define 

r^ 1 \p 1 ,p 2 ,p 1 ) = 2h{p) (21) 

By substituting Eq. ([21]) into Eq. ((20j) wc find 




Fig. 3 Diagrammatic identity for r^- 2 ' 1 ' in PA. 



Dividing i^ 2 ' 1 ) by r' 2 ' 2 ) one observes that their quotient is equal to A 3 /(2A4), 
so that the relation between three and four point vertices is not renormalized. 
Of course, for the actual PA model one must take A3 = 2A 4 = 2A, and 
conversely l(p) = h(p). In the rest of the manuscript, we only consider this 
case unless otherwise stated. 

In Appendix [E] it is shown, by using NPRG equations, that these ex- 
pressions are in fact non-perturbative (they are valid beyond an all-order 
perturbative analysis). These expressions have already been obtained before 
for the vertices with two incoming legs F5UI TU5T] . as a sum over bubbles. The 
interesting point is that the present analysis applies to any J^™'" 1 ) vertex in 
PA. As an example, in appendix |A"1 the equation for J^ 3 * 3 ) is obtained. Unfor- 
tunately, for n > 2 the corresponding equations must be solved numerically. 

Now that we have a method to calculate all correlation functions in PA, 
we can study BARW by means of a perturbative expansion in the branching 
rate a. We stress that a perturbative expansion on a coupling constant around 
a non-Gaussian model, such as PA, is a priori difficult to perform. 

To end this section, notice that all the results above are independent of 
the space dimension d. This allows us to make predictions independently of 
the upper critical dimension d c of the BARW systems studied below. 
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4 BARW - DP 

In this Section we consider the simplest BARW-DP model, where the only 
reactions are A — > 2 A and 2 A — > 0. More general cases in the DP universality 
class can be considered as well using the same methods. The microscopic 
action for this model reads, after the shift in the response fields (see Eq. ©) 

S DP = j (4>{d t - DV 2 )4> + X4>(4> + 2)^ 2 - a4>(4> + 1)0) (23) 

We now show how to perform a systematic expansion in a while keeping a 
finite A. This expansion is particularly well suited for properties of the model 
that take place at small a, but at values of A that can be out of reach of 
a perturbative expansion around the Gaussian theory. As mentioned in the 
introduction, the transition between the active and the absorbing phases in 
this model takes place, for d > 2, at values of A larger than a threshold X t h, 
which make the calculation of the phase diagram impossible within the usual 
perturbative analysis in these dimensions. As this threshold corresponds to 
a arbitrarily small, the value or Xth is computable in an exact way at the 
leading order of the expansion in a that we detail in the following. We stress, 
however, that the calculation of this quantity is just a specific example of 
an application of the expansion in er, which may be used for more general 
purposes. 

In order to analyze BARW-DP it is useful to consider, as in PA, a gener- 
alization of S DP with independent couplings. We then consider the action 



S DP = 



- DV 2 )cj) + A 3 00 2 + X A (U) 2 ~ ~ °30 2 4>) (24) 



As in the case of PA, one can deduce a Ward identity for the U(l) trans- 
formations Ijlip. which in this case reads 

ar dr r ( xr 7 sr\ n 

that leads us to 

(n - m)rM = X " a 3 ^- (26) 

aA.s do- 



where _r(™>™ 1 ) is a function of (x\, . . . , x n , x\, ■ ■ ■ , x m ). The solution of (j26|) 
implies the following relation for 7^( n ' m ) 

r ( "' m V2,<73,A 3 ,A4) = CT3 m -"7 ( "' m V2,<7 3 A3,A 4 ) (27) 

for m > n, with 

7 (n,m) a re 

gular function of its arguments (in particular 
for 03 = 0). This is nothing but the well known result of perturbation the- 
ory, which states, putting aside a re-scaling of vertices, that cubic couplings 
appear only via their product. At leading order in 173, Eq. (|2"T]) shows that 

r(«> m ) ~ G(af- n ) for n < m (28) 
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and that the calculation at leading order can be done at A3 = 0. 

In order to perform the 03-expansion one can consider the generating 
functional (|3]) and expand the exponential term 

Z = J V4>V4> exp ( - S DP \ CT3=0 + Jj^ + J^f^l ^ 3 jf 4> 2 <^J (29) 

In this way, the calculation to any order in 0-3 of any correlation function is 
reduced to the calculation of higher order correlation functions in a modified 
PA that includes a mass-like 02 term. It is worth mentioning that the methods 
presented in the previous Section work as well in the model including a 02 
term. When and if this 02 term is not necessary to make the theory IR safe 
it is possible to expand in 172 as well as 173 and this is what we are going to 
do in practice. 

As a final comment with respect to the cr-expansion, it is important to 
notice that it generates a convergent series, something not very common 
when dealing with perturbative expansions in field theories. This property 
follows from Lebesgue's dominated convergence theorem [55], given that we 
have under nonperturbative control the PA model (as shown in Appendices 
© O and [0). 

This convergence property can be most easily seen by working with a zero 
dimensional toy model 

Z = J dxe- Xx4+ax3 (30) 

Defining 

" 1 

/„(*)= ^Te-^-cr"* 3 ™ (31) 

i 

we see that the integrands 

< f n (x) < Ae~ Xxi (32) 

for some large enough constant A. From the dominated convergence theorem 
we know that the sequence of integrals of functions f n {x) converge to the 
integral of the function f(x) = liixin^oo f n (x). A similar reasoning applies in 
the case of the (d + l)-dimensional model, at least for the model defined on 
a lattice and in a finite volume. 



4.1 Threshold of the active-to-absorbing transition for BARW-DP 

Let us consider as a specific example the calculation of the threshold Xth 
for the existence of an active-to-absorbing phase transition in BARW-DP. 
Notice that this threshold value is non-universal, as would be a critical tem- 
perature in an equilibrium model. In order to check for the presence of a 
phase transition, it is enough to study the behavior of A — ^'^(p = 0) 
as a function of the annihilation rate A. In fact, we can detect this phase 
transition by looking for the zeroes of A, which correspond to a divergence 
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in the correlation length |30) . Given that X t h corresponds to the transition 
value of A when a — > + , an analysis at leading order in a allows for an exact 
calculation of Xth ■ Following the lines of the previous discussion, an equation 
for rW(p) at order 0(a) can be represented in the diagrammatic form of 
Fig. |4j that can be written 

r^\p) = -a + a f G(q)G(p-q)r^ 1 \q,p~q,~p)+0(a 2 ) 

Jq 

= -<r + Aal{p) I G(q)G(p - q) + 0(a 2 ) (33) 

Jq 

In the second line of (|33p we have evaluated the propagator G(p) and the 
vertex r ty2 ' 1 \q,p — q, —p) at order zero in a, and consequently replaced this 
last function by 4Z(p) (see Eq. (|19p . remember that we consider l(p) — U(p)). 




Fig. 4 Closed equation for F ( ' ' at first order in a in BARW-DP. 




Fig. 5 Closed equation for r (1,2) at first order in a in BARW-DP. 



As a side note, observe that we could have just as well written an equivalent 
equation for i^ 1 ' 2 ) (see Fig. [S]), which reads at order a 

r^ 2 \p 1 ,p 1 ,p 2 ) = -2a + 2a [ G(q)G( Pl ~ q)r^ 2 \ Pl - q,q,p 1 ,p 2 ) + 0(a 2 ) 

Jq 

= -2a + 8<r [ G(q)G( Pl - q)l( Pl ) + 0(a 2 ) (34) 

Jq 

As before, the r^ 2,2 ^{q,pi —q,pi,p2) vertex can be taken at order a , that is, 
it can be taken to be equal to 4i(pi). Expressions ([331) and (|3"4")) imply that 

r^ 2 \ Pl ,p 2 ,p 3 ) = 2r^\ Pl ) + 0(a 2 ) (35) 
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Which states that at first order in cr, the bare relation between the (1,1) and 
(1,2) vertices is maintained. 

Returning to our problem, we can look for a second order phase transition 
by studying the behavior of A. One needs the non-universal value l(p = 0) 
that can be obtained by evaluating Eq. (|l~9l) at p = 0: 



where 



^ = °) = TT^W (36) 

1(d) = f G(q)G(-q) (37) 

Jq 

By substituting the expression for l(p = 0), and evaluating (|33)l at p = one 
arrives at 

which for A = implies a threshold value 

A - = Wd) m 

To evaluate \ t h, we need to take into account that the properties of a phase 
diagram are not universal and depend on the specific form of the theory at 
small distances. This is as in equilibrium statistical mechanics, where critical 
temperatures depend on the specific form of the lattice. We will consider two 
particular microscopic forms for the model. The first one corresponds to the 
model defined on a hyper-cubic lattice with lattice spacing a. The second 
corresponds to a 'continuum' version where a UV cut-off is imposed at a 
finite (but large) scale A. 

For the hyper-cubic lattice, the propagator reads 

G(q) = 3-^ , (40) 

and the integral in ([55} becomes 



duJ f Ti^few-*) 

7T/a<qi<Tr/a l z7r ) 

d d q 1 



2 J- n/ a< gi <,/a W ^ Etl(l - COs(a 9i )) 



2-d 



1 



4£> (2vr) d ^ =1 (l-cos( % )) 



(41) 



where the integral over w has been performed by using the residues' theorem. 
The remaining integral must be calculated numerically. In Table [1] the value 
of the resulting threshold coupling is given. Previous results from Monte- 
Carlo simulations and approximated NPRG equations jl6][3T1 are in good 
agreement with these exact ones. This same general structure of the phase 
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d 




3 4 


5 


6 


Xth/Da d - 
\th/Da d - 


~ 2 (this work) 

" 2 (Monte-Carlo) [16] 


3.96 6.45 
3.99 6.48 


8.65 
8.6 


10.7 
10.8 



Table 1 Values of the threshold coupling Xth for various dimensions d. Comparison 
of present exact results with Monte-Carlo [16] . 



diagram has also been shown to exist in other models in the DP universality 
class [52] . 

An interesting property observed in [TB] is that Xth seems to grow linearly 
with d. In |33| . a single-site approximation scheme that is argued to become 
exact in the large d limit in a hyper-cubic lattice was analyzed, and this 
linear behavior was obtained. In order to analyze such a behavior here, it is 
necessary to find the large-d limit for the integral 1(d). For this purpose it is 
useful to re-write it in the following form: 

I(d)Da d - 2 = — [ ^ (42) 

To solve it, one can imagine the various cos(gi) as random variables with 

zero mean. By the strong law of large numbers, their mean ( J^i=i C0S ('Zi)) /d 
tends to zero, except in a zero measure set. We are then tempted to substitute 
the limit inside the integrand and obtain 

I(d)Da d ~ 2 d ~°° ^. (43) 

This step is non trivial from a rigorous mathematical point of view, but 
turns out to be correct by using elaborate methods of real analysis Q- As a 
consequence, 

X th /Da 2 - dd ^°2d (44) 
in agreement with previous results [33) . 

It is interesting to observe that expression (]39[) only depends on quantities 
that are calculated exactly in the Local Potential Approximation (LPA) of the 
NPRG, which is the lowest order of the Derivative Expansion. Only vertices 
at zero momenta are used and their exact equation turns out to be the same 

1 In order to prove this property in a mathematically rigorous way, one needs to 
prove that if z n = l/ n X) n cosg n ,then f dq\ . . . dq„(l/(l — z„))° (*) is bounded by a 
constant independent of n, with a > 1. This would show that the sequence 1/(1 — z n ) 
is uniformly integrable and that f dqi . . . dq n l/(l — z n ) —¥ 1. In order to find a 
bound for (*) one can use that f dq\ . . . dq n (l/(l — z n )) a = J °° dxProb((l/(l — 
z„))° > x) < xo + dxPiob(z n > 1 — x _1//a ) (**) where the probability is taken 
in the sense of considering the q n as random variables uniformly distributed in 
the interval [— ir,ir]. In order to give a bound of the second integral in (**), one 
can use standard techniques of large deviations bounds in order to prove that 
Prob(jz n > 1 — x~ 1/,a ) < l/x nC / a with a constant C independent of n. This gives 
a bounded integral for n> a/C, proving the desired result. 
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as the one that comes from the LPA (see Appendix |B|) . This a posteriori 
explains the success of the LPA in reconstructing the phase diagram of this 
model [IB] . 

However, as mentioned before, the phase diagram is a non-universal prop- 
erty that depends on the precise definition of the model in the ultraviolet. 
In particular, the value of the integral 1(d) is different if calculated in a 
discrete lattice or in the continuum with a given ultraviolet regularization. 
In the previous study done within the NPRG [16] . a continuum version of 
the model was implemented, but the initial bare condition was imposed at 
a finite (but large) value of the microscopic scale A, which serves as a UV 
cut-off. In order to be able to compare our continuum non-universal results 
with those obtained in [16] we will choose an UV regularization compatible 
with the NPRG procedure, as described in Appendix iBl 

In the continuum regularized case, the integral to be calculated in order 
to make a direct comparison with the NPRG (see Appendix [Ej) is (with the 
tilde indicating this second regularization) 



1(d) ' 



(4n) d / 2 r(d/2) D 

A d-2 



(45) 



(4ir) d / 2 r(d/2) D d(d — 2) ' 
This yields for this particular regularization 

- A 2 - d (4n) d / 2 r(d/2)Dd(d-2) 

x th = ^ (4b) 

Given that this integral is calculated in a closed form by using exclusively 
quantities evaluated at momentum p = 0, we can check (see Appendix|E]) that 
it coincides exactly with the LPA equation for this same quantity. Indeed, 
our result recovers the numerical LPA solution of the NPRG of [TB] within a 
nine digit accuracy. 

We can also compare the results coming from both lattice and continuum 
regularizations, as has been done in |161I34) . As is explained there, one cannot 
do such comparison without fixing the relation between A and the lattice 
spacing a. In [TB] , this relation was fixed by multiplying the continuum results 
by exp(c(2 — d)) and fitting the constant number c, finding a very reasonable 
agreement up to d = 7. However, we observe in the present results that the 
agreement is lost in higher dimensions where the continuum version leads to 

-^^^Me^^ (47) 

This indicates that the agreement between both results is only valid for a 
limited range of dimensions. In order to relate the results in a larger range of 
dimensions, one must consider a d-dependent relation between A and a or, as 
done here, take into account the precise ultraviolet regularization considered. 

Finally, it is convenient to point out that for d < 2 an IR divergence of 
the integral in (|39[) takes place. This makes Xth — in those dimensions, in 
agreement with the results of [8]. For this reason, for d < 2 it is not useful 
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to expand the model at small a for a finite A in order to study the phase 
transition. Moreover, this also shows that in those dimensions the transition 
is dominated by IR effects, and correspondingly most of the dependence on 
the microscopic behaviour of the model is absent. 



5 BARW - PC 

Let us now consider BARW-PC, corresponding to the Parity Conserving/Genteralizcd 
Voter universality class. In this case, it is convenient not to shift the response 
field in order to make explicit the (f> — > —(f), (/) — » — <j) symmetry associated 
with conservation of the parity of the number of particles. The microscopic 
action for the BARW-PC model reads (see Eq. ^) 

S PC [^, $]= [ U{dt - DV 2 )cf> + A(0 2 - l)<f> 2 + <r(l - 2 )#) (48) 

J X 

where the last term corresponds to the branching reaction A — > ZA with rate 
a. 

The case a = corresponds to Pure Annihilation, now written in terms 
of the non-shifted </> field. This version of Pure Annihilation can again be 
solved following the same ideas as previously. Here, as opposed to the shifted 
case, we have the additional constraint that r^ n ' m ' = if (n + m) is odd. 

Let us now show that in this version of PA 

r("' m ) ~ ( A ("-™)/2) forn>m (49) 

and zero otherwise. We again define a generalized action S with indepen- 
dent A2 and A4 couplings as in Eq. (ITUl) 

S PC [(j )l 4>]= I U(d t -DV 2 )0-\ 2 cj ) 2 + \ 4 $ 2 ( j ) 2 + cT 2 4, ( j ) -a 4 4> 3 A (50) 

First we set 02 and 03 equal to zero, in order to be in PA, and exploit the 
U(l) Ward identity for the infinitesimal transformation 

4>(x) -> (l + e)4>(x) 

$(x)->(l-e)fa) (51) 
The argument is completely analogous to the one shown in Section[3j yielding 

(n - m)r^ = 2A 2 ^— (52) 

(with, as before, J^"*" 1 ) a function of (x%, ■ ■ ■ , x n , x%, . . . , x m )) from which 
Eq. g9]) follows. 

It is easy to check that the equation for J^ 2 * 2 ) remains the same as in the 
shifted case, Eq. (|17l) . and we thus define the function l{p) again by means 
of Eq. (1191) . The vertex 

r (2,o) 

can be studied by following similar lines, and 
is found to be related to l(p), by r (2 ^(p) = -2l(p). Also as before, ft 1 ' 1 ) is 
easily proven not to be renormalized in this version of PA. 
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Since we are interested in studying the cr-expansion around PA, it is useful 
to establish the equivalent of Eq. (051) regarding the order in a of the r^ m \ 
We again work with generalized couplings a 2 and 04, using the modified 
action (|50[) and we arrive, by using the U(l) Ward identity deduced from 
([5Tjl (an identity similar to Eq. ([25])) at the relationship 

r^ m \a 2 ,a 4 , A 2 ,A4) - ^""^V"'^^,^, A 4 ) (53) 

for to > n, with gular function of its arguments (in particular for 

cr 4 = 0). This implies that r("> ro ) - 0(a^ n -' n ^ 2 ) if to > n. The details of 
the calculations leading to this property are completely analogous to those 
of Section H in BARW-DP. 

5.1 The stability of the PA fixed point 

One striking feature of the PC model is the existence of an active-to-absorbing 
phase transition in d = 1 . This is believed to be related to a change of stability 
of the PA fixed point in a dimension d c between one and two. Perturbatively, 
and also within the LPA, this change of stability occurs in the following way 
(see a schematic representation of this scenario in Fig.lo1[8.ll9j). On one hand, 
in d = 2, the Gaussian and PA fixed points merge so that, for dimensions 
close to two, the relevance of the branching reaction A — > 3A can be proven 
by canonical power counting arguments. On the other hand, at 1- and 2-loop 
orders an (upper) critical dimension d c > 1 is found such that for d < d c 
the coupling a becomes irrelevant around the PA fixed point which there- 
fore becomes fully attractive. This change of stability occurs because a new 
fixed point, F PC , crosses the PA fixed point at d c and in this dimension they 
both change their stability. Below d c , this new fixed point is in the physically 
relevant quadrant A > 0, a > 0, has one unstable direction, and is thus asso- 
ciated with the phase transition. The PA fixed point is then fully attractive 
for d < d c and describes the absorbing phase. Notice that the value of d c 
changes significantly between 1-loop - where d c = 4/3 - and 2-loops where 
d c ~ 1.1 (d c = 4/3 within the LPA). 

Some of these facts seem to be confirmed by other methods. In d = 1 
Monte-Carlo simulations of this model show indeed a new universality class 
38,37 , and an exactly solvable model expected to be in the same universality 
class as BARW-PC shows a negative scaling dimension for a: d a = —1 35 . 
This result (d a = —1) is identical to the prediction at order e [5]. At two 
loop order, though, this value for d a changes and gets smaller in magnitude, 
d a ~ — 0.137 at d= 1 [5]. As will be argued in the following, this significative 
difference between MC and 2-loop results can be seen as a first indication 
that the results of [35] are not entirely valid for this system. 

We now reanalyze the stability of the PA fixed point in the presence of 
the PC creation reaction, A 3A, that we can determine exactly since our 
analysis is exact at small a. The relevance of this coupling can be obtained 
from the flow of either J^ 1 ' 1 ) or r^' 3 \ since both these functions are of order 
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a 



A active phase 



absorbing pjhase 



d, 



2 d 



Fig. 6 Sketch of the relevance of th e a perturbation in BARW-PC around the PA 
fixed point, as expected from 8,19]. The arrows show the direction of the RG flow 
for the coupling a. Above d c , a is relevant, whereas it is irrelevant below d c . The 
dashed line represents the location of the fixed point F PC that crosses the PA fixed 
point at d c and that is associated with a phase transition below d c . 



a. However, the RG flow of r^ 1 ' 1 ' depends on the somewhat difficult to study 
p(3,i) ver t ex f pa ( se e Fig. [7]), and we prefer to study r^ 1 ^ . 



+ 





Fig. 7 Closed equation for r (1,1) in BARW-PC, at first order in a. 



At first order in a, any diagram for J^ 1 * 3 ' will be of the form shown in 
Fig. |SJ As can be seen in this figure, it involves the bare a vertex as well as 
the PA r( 2 ^ and I^ 3 ' 3 ' 1PI vertices. As it stands, though, we would have 
to solve the independent equation for _r( 3 ' 3 ) in order to make progress (an 
analysis of which can be found in Appendix [21 where it is shown that its 
equation requires numerical methods to be solved). Moreover, this expression 
is not well suited for the analysis of universal properties, because it is ex- 
pressed in terms of the bare vertex and not in terms of the full i^ 1 ' 3 ) vertex. 
Fortunately, the fact that we only deal with PA vertices allows us to find an 
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Fig. 8 Closed equation for r (1 ' 3) in BARW-PC at first order in a. 



easier relationship for r^' 3 \ using the already known property which allows 
to find closed forms for PA vertices. Notice that J^ 2 - 2 ) and I^ 3 ' 3 ) have always 
two possible closed decompositions, being vertex of the form J^™*™) (see dis- 
cussion in Section [3]). Specifically, we can rewrite the equation for J^ 1 ' 3 ) in 
the form shown diagrammatically in Fig. [9l which can be written 

- 2A f G(q)G(p 1 + p 2 - q)rW (p, q, Vx + P2 - 9, Pa) 

Jq 

- 2A f G{q)G{p x + p 3 ~ q)r^ (p, q, p x + p 3 ~ q, p 2 ) 

Jq 

-2A [G(q)G(p 2 +p 3 -q)r^(p,q,p 2 +p 3 -q,p 1 ) (54) 

Jq 




The highly symmetric form of this equation suggests the following ansatz for 
the functional form of Ja 1 ' 3 ) (which can be easily checked using the equality 
p=Pi +P2+P3) 



r (1 ^ ) (p,p 1 ,p 2 ,p 3 ) = -2<r(p,pi) - 2<r(p,p 2 ) - 2<r(p,p 3 ) 



(55) 
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In terms of cr(p,p) the equation becomes 
a(p,p)=a-2X [ G(q)G(p-p-q) 



x (o-O, q) + a(p, p - p - q) + <r{p, p)J (56) 
Using the known expression for l(p), Eq. (flU]) . we obtain 



<j(p,p) = jl(p-p)-2l(p-p) I GKDGin- j, ,,) 



x fo"(p, <7) + <x(p,p - p - q)J (57) 



For the calculation of d c , it is enough to analyze the p — behavior (we are 
interested in the IR fixed point structure of the theory). Defining 

trip) = trip = 0, -p) (58) 

(notice that we have chosen a minus sign in the definition), and after a change 
of variables inside the integral, we get 

trip) = ~J(p) ~ 4J(p) ^ G(q)G(p - q)u{q) (59) 

From now on we will omit the tilde in p for notational simplicity. The quantity 
we are mostly interested in is d a , the scaling dimension of a in the IR limit 

(rip) - \p\ d " for u, |p| 2 < A^J (60) 

Naive power counting yields d a = 2, which would imply that the branching 
perturbation a is relevant for all d, but fluctuations of course change this value 
of d a , and could even make it negative, which would imply the irrelevance of 
the a perturbation. 

In order to solve Eq. f|59[) it is useful to define the quantity 

<»> 

whose behavior in the IR is expected to be of the form a(p) ~ |p| d_£iCT (recall 
that l(p) ~ j j-, j 2 rf - n fag^ re gi me )_ The equation for a reads 

Hp) = j- A f G ii) G (p - i)Hi)i(u) (62) 



Using this exact expression and expanding in e — 2 — d we recover the 1-loop 
result d c — 4/3, as well as the 2-loop result d c ~ 1.1 [5]. These results follow 
from a perturbative series in a and A, and from a simultaneous expansion in 
e = 2 — d- The details of these calculations can be found in appendix |FJ 

In order to get an exact result for d a it is convenient to get rid of the bare 
reaction rates, as we are interested in the universal IR scaling behavior. Let 
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us start by doing so in the case of l(p), which will be useful in what follows. 
The IR limit is taken by making A — ¥ oo (more precisely, by considering 
v i Ip| 2 ^ \ 2 /( 2 ~ d \ the typical momentum scale set by the bare annihilation 
rate). This can be done safely for d < 2, and is a subtle limit when one 
studies directly d = 2 in order to perform the e-expansion, as will be further 
commented in Appendix |FJ By exploiting expression (|19[) we obtain the IR 
behavior 



{l IR (p))- 1 = 2 [ G(q)G(p-q) 

Jq 



J (2ir) d J 2ir q 2 + iuj (p - q) 2 + i(y - u) 1 } 

and thus 

Now we can return to a . As we are only interested in its scaling behavior, 
it proves convenient to subtract to (|62p its value at zero a(p = 0), which is 
zero in the IR for d < 2, given that we expect d„ < d. This is seen to be true 
in the e expansion around d — 2, and must be true near the sought-for d c , 
where d a should be zero. Our results will later confirm d a < d. We thus have 

HP) = - 4 J q °(lW)G(q) (g( P -q)- G(-q)) (65) 

This is a complicated equation, and to be able to solve it, we must take into 
account the scaling invariance we expect from its solution. We exploit scale 
invariance in order to define the scaling function d(y) 

a(p,v) = \p\ d - d ^(v), u=^ (66) 

Observe that we are performing a perturbation around the PA fixed point, 
whose anomalous dimensions are zero (that is, 77 = 0, z = 2, as already 
mentioned). Accordingly, the natural scaling variable is v = i//|p| 2 . 

We can now write an equation for a(C>), using the form (|64p for l(p) and 

choosing as variables Q — cu/q 2 , q = \q\/\p\ and u — cos (p, q) 



a{v) = -4 



d — 1 



(4^) rf / 2 \ / 2tt^ 



2 1 -<V 2 r(i - d/2) J \{2Tr) d + l r{^ 

I dqq d - d " +1 [ du(l-u 2 )^ I <IJj a(Sc) [ - + /J: 



1 f oo / j \ l-d/2 



n J -1 J-00 \2 

1/1 1 



1 + iui \ 1 + iv + q 2 (1 + iijj) — 2qu q 2 (1 + iuj) / 



(67) 
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Re[a] 




Fig. 10 ^-dependence of the real part of the scaling function a(y) for several values 
of d. 



This equation is still too complicated to be solved analytically, and requires 
a numerical solution. A convenient way to do that is to make an expansion 
in u, which turns out to be rapidly convergent. We then proceed as follows: 
at each order in the expansion in u we adjust da- at a given value of d, by 
numerically iterating this equation in order to reach a fixed functional form 
for a(y) in a lattice of N v points with a resolution 5v. We have checked the 
convergence in u and in the numerical parameters 8v and N„, used for the 
computation of integrals. This procedure gives always a converged scaling 
function cr(C>), which confirms a posteriori the scaling form ansatz (I66[) . In 
Fig. [TU1 we show the explicit v dependence of the function 5{v) for some 
values of d. As can be seen, it is a non-trivial function of its argument, 
which may explain the qualitative difference between our results and previous 
approximate results. Observe that LPA and 1-loop analysis are based on a 
constant coupling a (without dependence on frequency and momentum). As 
expected, this dependence becomes weaker as d approaches 2. 



This procedure allows us to find the value of d a as a function of d, the 
results of which are plotted, together with previous perturbative results, in 
Fig. [IT] There one can see that even if d a gets smaller when d decreases, 
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d 



Fig. 11 Results for d CT , showing there is no change in the RG relevance for the 
branching rate a. 



it remains always positive. This is an unexpected result, which deserves a 
careful discussion. First of all, it is important to observe that this result does 
not rule out the existence of a new fixed point F PC for small d, corresponding 
to the PC universality class. A new fixed point can indeed appear but for a 
nonzero value of the branching rate a*, as seen for example in the sketched 
flow shown in Fig. Q21 This would mean in particular that the low branching 
phase of the model would have a behavior different from PA. This can be 
studied either by using Monte-Carlo methods or perhaps by means of the 
NPRG at orders higher than the LPA. 

As said before, there exists an exact result in d = 1 [35] which seems to 
indicate that in this universality class a is in fact irrelevant with respect to 
the PA fixed point in that dimension. We can explain this difference observing 
that the model used in [35] is defined with A = oo, and indeed presents no 
phase transition at all for whatever value of a. Now, the IR limit corresponds 
to v, |p| 2 <C \ 2 /( 2 ~ d \ but this does not allow us to take a = when compared 
to A. Looking at Eq. (p2")l . A = oo implies a = 0, so that the relevant direction 
corresponding to a is no longer accessible by studying a as a perturbation. 
This is true for all d. Indeed, the results of [5] are also compatible with this 
scenario: in fact they show an irrelevant a for all d when A = oo. For example, 
at 1-loop level (approximation valid close to d — 2) we have (see appendix 
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Fig. 12 Sketch of the relevance of the a perturbation in BARW-PC, compatible 
with the results in this work and with simulations. 

ED 

Hp) - ^i 2 (p) (68) 

so that we see explicitly that A — > oo yields a = 0, and a is no longer associ- 
ated with the relevant branching direction. This is in evident contradiction 
with the fact that there is an unstable direction in that dimension. Thus, we 
think that the exact calculation in [55] does not apply to BARW-PC, the 
system in which we are interested. 

Also, there exist a result in |36) in which branching and annihilating 
systems of particles performing Levy flight dynamics are studied. In it, the 
authors show that a change in the value a of the Levy flight exponent can 
be made to correspond to a change in the dimension d of the corresponding 
standard BARW system. This is used to recover d c = 4/3 for BARW-PC. The 
analysis, however, is made by means of a 1-loop perturbative expansion in A 
and <7, which explains the coincidence with the results of 0. A re-analisys 
of Levy flight dynamics can be envisaged within the approach propposed in 
this work. 

In what respects Montecarlo studies of the low branching regime of this 
system, they have until now, as far as we know, also been mostly made in 
the limit A — > oo [571155] mentioned before. They are compatible with the 
standard scenario, but within the criticims previously pointed out. 

6 Summary and discussion 

In this work we have applied field theoretical methods to answer some non- 
trivial questions about a class of reaction-diffusion systems. We have pro- 
ceeded by exploiting the special case of Pure Annihilation, a system which 
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does not present a phase transition but which nonetheless possesses a non- 
trivial fixed point in the RG sense. In order to do so, we took advantage of 
its simple structure, as well as the symmetries and causal properties of the 
system (which in fact allowed us to go beyond perturbation theory) . We have 
then applied an expansion in the branching rate around Pure Annihilation, 
giving us access to the small branching regime of BARW, both with and 
without an additional parity conserving symmetry. 

We have chosen to concentrate, as a first order example, in some im- 
portant properties of these systems, usually very difficult to control but that 
become possible to solve within the present method. In the case of the system 
of reactions 2A — > 0, A — > 2A, which belongs to the DP universality class, we 
have given an explicit proof of the existence of a phase transition in all space 
dimensions, already seen in previous numerical solutions of approximated 
versions of the NPRG flow equations, and in Monte-Carlo simulations. We 
have moreover calculated exactly the non-universal threshold value for the 
annihilation rate in order to find this phase transition in two sample systems. 
This result is beyond the possibilities of usual perturbation theory. 

In BARW-PC, where the parity of the number of particles is conserved, 
we have concentrated on the value d c of the upper critical dimension, that 
was previously believed to be somewhere between d = 1 and d = 2. Previous 
1-loop and LPA results indicated d c ~ 4/3. By truncating our equations at 
one- loop order we were able to recover this approximate result, as well as 
the two-loop result of [5]. Surprisingly, we have found that the appearance 
of the PC fixed point associated with d c must occur at a nonzero value of 
the branching rate, which would be compatible with a scenario where there 
exists not one but two new fixed points for d < d c . Further investigation of 
this issue should be performed, either by a higher order expansion in a or by 
lattice simulations, or by the use of the NPRG method at orders higher than 
the LPA. Work in some of these directions is already underway. 

Let us emphasize that the cr-expansion here introduced represents an 
expansion around a non-trivial (non- Gaussian) model, which in particular 
implies, as explained in the text, that the first order results obtained in this 
work represent the first term in a convergent series. 

The cr-expansion has allowed us to obtain results not accessible with the 
usual perturbative expansion, while still being (for the most part) analytical. 
This kind of ideas could in principle be generalized to other field theoretical 
systems, and future applications can be thought within the study of out of 
equilibrium systems. In what concerns BARW, a second order expansion in a 
would in principle allow for the approximate calculation of critical exponents. 
Extensions to other out of equilibrium systems can also be envisaged, for 
example in the study of PCPD, or of the Cole-Hopf version of the KPZ 
equation [27] . which has a structure very reminiscent of Pure Annihilation. 
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A Study of the J^ 3 ' 3 ) vertex in PA 

In this appendix we will study the _r( 3 - 3 ' vertex in PA with the methods intruduced 
in Section [3] Unfortunately, as we will see, this will not be enough to find a com- 
plete analytical solution for the vertex. Such a solution would need a numerical 
implementation, beyond the scope of the present work. 




+perm. 



Fig. 13 Closed equation for r (3,3) in PA. 



As explained in the main text, every diagram for _r' 3 ' 3 ) is of the form shown 
in Fig. [T] If we analyze the possible 1PI contributions we end up with the two 
diagrams shown in Fig. [T3] Using the known form for l(p) and the corresponding 
symmetry factors, the equation corresponding to this diagram reads 

r (3,3) (pi,P2,P3,Pi,P2,P3) = 

64A ^ I G(q)G( Pl +p 2 -q) [l(p2 + p:i)G(p 3 + q- pi)l(q + p 3 ) 
+ l(pi + Ps)G{p3 +q- P2)l{q + Pa) +l{pi+ P2)G{p 3 + q- pa)l(q + pa) 
+ f G(q)G( P i+ P 3-q)l(p 2 + q)\G{ P 2 + q-pi)l(p2+P3) 

Jq 

+ G(p 2 + q - P2)l(pi + Ps) + G(p2 + q- P 3)l(pi + pa)] 
+ / G(q)G(p2+P3-q)l(pi+q)\G(p 1 +q-p 1 )l(p2+p 3 ) 

Jq 

+ G(pi + q - p 2 )l(pi + P z) + G(pi + q — P-a)Kpi + P2)] 

- 2A^G(g) (g( Pi +P2- q)r (3 ' 3 \q, Pl + P 2 - q, P3 ,pi,P2,p 3 ) 
+ G(pi +p 3 - q)r (3 ' 3) (q,p 1 +P3- q,p2,pi,P2,Ps) 

+ G(p 2 +P3 - q)r (3 ' 3) (q,p 2 +P3 - q,pi,pi,p2,p 3 ) ) (69) 



This equation shows a very symmetric structure, which suggests the ansatz (easily 
proven to be correct using momentum conservation) 

r (3 ' 3 \p 1 ,p2,P3,Pi,P2,P3) = f(pi +P2,P3,Pi +P2,P3) + permutations (70) 
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where in fact the last dependence (in p3 in the equation) is redundant due to 
momentum conservation. The corresponding equation for / is 

f(Pa,Pb,p a ,Pb) = 64\l(p a ) / G(q)G(p a - q)G(p b + q - p~b)l(pb + q) 

Jq 

- 2A / G(q)G(p a - q)\f(Pa,Pb,Pa,Pb) + f(q+Pb,Pa - q,Pa,Pb) 

Jq 

+ f{Pa +Pb - q,q,Pa,Pb))\ (71) 
and, using the explicit form for l(p) this can be rewritten as 

f{Pa,Pb,Pa,Pb) = 64:l( Pa )l(p a ) / G{q - p b )G(q - pb)G(pa + Pb + q)l(q) 

J q 

-4l(p a ) / G(q)G(p a - q)f(p a + Pb - q,q,Pa,Pb) (72) 



where a change of variables has also been performed inside the integrals. 

In order to get rid of the incoming and outgoing dressed l(p) vertices we can 
define 

, ~ ~ s fiPa,Pb,Pa,Pb) ,„„s 

X{Pa,Pb,Pa,Pb) = ~ - (73) 

l{Pa)l{Pa) 

whose equation reads 

X(Pa,Pb,Pa,Pb) = 64 / G{q - p b )G{q - Pb)G{p a + Pb + q)l{q) 

J q 

I / G{q)G{p a - q)l{pa + Pb - q)x(Pa + Pb - q,q,p a ,Pb) (74) 

J Q 
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This expression can be further simplified, but in the long run a numerical study 
of it is unavoidable in order to solve it. Notice also that this expression is not 
explicitly symmetric with respect to the change p — ^ p. It is easy though to find 
such a symmetric expression, by combining this expression with the one obtained 
by using the equation stemming from the diagrammatic ansatz with the A bare 
vertex in the outgoing legs of the diagram. This symmetric equation does not turn 
out to be simpler to solve than (f74)l . 



B Non Perturbative Renormalization Group study of BARW 

The Non-Perturbative Renormalization Group [31?ll4TJir4TTl42U43ll4"4"] is a general frame- 
work for the study of strongly correlated systems, mainly for field theoretical prob- 
lems. It is basically the modern version of Wilson's initial ideas on the RG. One of 
the main interests for using the the NPRG is that it enables to devise new approx- 
imation schemes, useful to tackle problems where more classical tools are of little 
or no use. In particular, the NPRG stands out as a natural ground for the study 
of strongly coupled systems, and the last years have seen an increasing number 
of applications of NPRG methods to d iver se kinds o f problems in many areas of 
(in general str ongl y corr elated) physics [39 45 , 46 , 47 , including out of equilibrium 
systems 4,5,6,16,1948 . In this regard, the NPRG is particularly well suited for 
the study of critical scale-free regimes, where renormalization group methods are 
known to give a simple description of the critical behaviour, and where strong cor- 
relations are gen erally p resent. Detailed introductions to the NPRG can be found 
in the references 39 , 30 , and also in particular when applied to out of equilibrium 



problems in [411 341146] 
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The NPRG formalism relies on the construction of a sequence of scale-dependent 
effective models for the model under study, interpolating between the short scale 
physics at the assumed microscopic scale A and the full long distance physics when 
the running scale goes to zero, by averaging over fluctuations at each value of this 
sliding momentum scale k. Instead of working at the level of an effective Hamil- 
tonian, as in the original Wilsonian approach, the NPRG is developed in terms of 
the flow of effective average actions which are scale dependent modifications 
of the 1PI generating functional r, in such a way that only takes into account 
fluctuations with characteristic momenta \q\ > k. At the scale k — A, no fluctuation 
is taken into account and Fa coincides with the microscopic action S, whereas full 
system information is recovered in the limit /&_►() = F. The procedure for con- 
structing the effective average action Gammak consists in the addition of a scale 
dependent mass- like term which partly freezes out the slow modes. This is achieved 
by adding a regulator term to the original microscopic action of the form [39] 

AS h [<t>A]= j (<fi(x),$(x))-R k (V 2 ,d t )-(<j>(x),$(x)). (75) 

In the expression ()75|) . Rk is the (model dependent) regulator matrix, with matrix 
elements ~ Rk, a cutoff function which behaves as Rk(a^,u]) ~ k 2 (in Fourier 
space) for small momenta |q| < k, in order to generate a decoupling of the "slow" 
modes, and Rk vanishing for large momenta |q| > k, so that the rapid modes 
remain almost unaltered. The scale-dependent generating functionals Zk [J, J] = 
J T>cj>D(j> exp( — S — ASk + J x J(j> + f J 4>) axe then used to obtain the effective 
average action Fk , through the Legendre transform of log Zk [J, J] ■ 

At the core of the met hod there is an exact functional differential equation, 
the Wetterich equation |44| . which governs the RG flow of Fk as a function of the 
scale k. This equation can be cast in a way that makes explicit its one-loop, 1PI 
structure: 

^ = 5^/0^^^) (76) 

where dk is the derivative with respect to k applied only to the explicit dependence 

" (2) 

on k of the regulator profile Rk, and F^ is the matrix of second derivatives of 
_Tfc with respect to the fields. One observes that the right hand side of Eq. (|76p 
is similar to the 1-loop expression for the r functional, but with vertices and 
propagators extracted from r k + ASk instead of S. We make extensive use of this 
1-loop structure in the following. 

Equation (|76[) cannot in general be solved exactly, and one usually has to per- 
form some approximation in order to extract physical inform ation fr om it. The most 
broadly used method is the so-called Derivative Expansion [411143] . which imposes 
a particular ansatz for the functional form of in the small momentum regime 
(in fact, for critic al syst ems it can be argued that it is valid only in the limit of 
momenta |q| < k 49 50 ). The DE, as well as other possible approximations used 
in the literature 39 , 49 , 50 , does not rely on any small coupling parameter, so that 
the approach remains non-perturbative in essence. It is of course also possible to 
perform usual perturbative calculations within the NPRG formalism. 

Here we assume the existence of a suitable regulator function, and in particular 
one that does not break the symmetries of the field theory. There is still a great 
degree of arbitraireness in the selection of such a regulator. 

As stated in the m ain text, causal properties, and in particular the so-called 
Ito prescription 25, 26 are much less trivial to study in the NPRG formalism when 
compared with the usual perturbative approach [3] . One can show that the earliest 
time in the time dependence of a vertex must correspond to a <j> field, and conversely 
the last time must correspond to a ^ field. This implies that the 2-point function 
Jfc remains causal for all k. The proof of this causal property will be shown in 
appendix [D] 
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In the particular case of PA, it will be shown in Appendices [C] and [D] that this 
property is generalized to 3 first-in-time <f> legs and 3 last-in-time (f> legs for any 
rj_ n,n \ This implies in particular that the 4-point function rj, 2 ' 2 * 1 and the 6-point 

(3 3) 

function ' are completely causal, in the sense that all its incoming legs must 
correspond to times smaller than those associated with any of its outgoing legs. As 
we shall see below, this property greatly simplifies the study of vertices with a low 
number of legs. 



C Non-perturbative proof of the form of the closed equation for 

r (n,m) j n p A 

In this Appendix we will use the NPRG as discussed in Appendix [B] in order to 
show the general diagrammatic property used in the main text (see Fig. [TJ , which 

relates i^( n > m > with A and vertices of at most (n, m) legs. This relation has been 
presented in the main text at all orders in a perturbative expansion in A, but here 
we sketch a non-perturbative proof for it. 

First, observe that we can write the property we want to show in terms of a 
specific form for F k , the average effective action for the theory at any scale fc 



r k = s PA + 



( <t> 2 {x) f G k (x'-x)G k (x"-x)r h [x',x",<f>,4>] (77) 

J x J x' ,x" 



with 



<5r fc [x , ,a: ,, ,^,< 
8cj>(x) 

sr k [x',x",cp, t 

5(j>{x) 



if t < max(t',t") 

if* < max(t\t") (78) 



indeed, this form ensures that any r™' m vertex will consit in a series of bare terms 
given by the action S PA , together with renormalized terms which always begin 
(temporally speaking) by a bare (2, 2) vertex. 

The property we want to show is obviously true at the bare level, with Fa = 0. 
Now we proceed in an iterative way, assuming that the property we want to prove is 
true at a RG scale ko, and checking that it continues to be valid for a scale ko — 8k. 

By hypothesis then, we take that at k = fco all vertices r^ n can be decom- 
posed as a diagram containing a A4 bare vertex (with two simultaneous incoming 
<j> legs) and a r^ n ' m \ this last function being constrained by the condition that the 
full diagram must be 1PI. As explained before, the NPRG equations for any ver- 
tex r^"' m) can be represented diagramatically by 1-loop diagrams, where vertices 
and propagators are read from r k + AS k - To those diagrams one must apply the 
operator d k in order to obtain c^r^™'™'. Each one of these terms will consist in 

a number of rj, 1 ' vertices joined together by an internal loop of propagators. We 
can distinghish between internal lines, pertaining to the internal loop, and external 
lines. 

Consider then a typical diagram contributing to d k r^ l ' m ' > (ti, . . . , t n , ii, • • • , t m ), 
where we emphasize the time dependence of the vertices. We now define as to the 
smallest time for any incoming external leg, to < tt Vfc. Its corresponding leg is 
attached to one of the vertices in the loop, and by hypothesis at k = fco, the two 
smallest times in this vertex must correspond to an incoming bare A4 vertex. We 
have then two possibilities to consider: 

1) Both to legs are external legs. This means that d s r^ n ' m ' 1 can also be decom- 
posed in the form ()78|l . Which implies the desired property for fc = ko — 8k 
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2) Only one to leg is an external legs, and the other one is an internal incoming 
leg. Thus, given the non-renormalization of the propagators, there is another vertex 
with an outgoing leg at a time previous than to- But at k — ko this vertex should 
have at least two incoming legs with a corresponding time previous to to- At least 
one of these incoming legs must be external, contradicting our assumption. 

The desired property is then preserved all along the NPRG flow, showing that 
it is a fully nonperturbative result as anounced in the main text. 



D Causal properties in BARW and PA within the NPRG 

In this appendix we will describe the non-perturbative causal properties in BARW 
and PA. Altough these properties have a simple expression in terms of a perturba- 
tive expansion [HJ, they are somewhat difficult to prove in the context of the NPRG 
[4] , which is what we choose to do here, in order to have non perturbative versions 
of the results. 

The method used here is completely analogous to that used in appendix [C] 
in order to prove the diagrammatic property used in the main text. In fact, this 
property can be seen as a causal property of PA. Here we show that it can be 
slightly generalized. But before that let us start with a simple causal property for 
generic BARW. 




Fig. 14 An example of a diagram contributing to 9fe_T^ n ' m in the NPRG for a 
generic theory. 



Let's concentrate on theories where at bare level, all vertices verify that in the 
time domain 

rjC' m) (t 1 ,...,t n ,i 1 ,...,i m ) = o (79) 

if one of the times in the set {ti, . . . ,t n } is strictly larger than all times in the set 
{ti, . . . ,t m }. For the 2 point functions this is equivalent to imposing the Ito pre- 
scription 25,4 by requiring that the propagator in time is strictly causal (G(t,t' = 
i) = 0). For all other correlation functions this can be seen as a consequence of the 
locality of the bare theory. These properties apply for all BARW, as can be checked 
by looking at Eq. (J6]l . We are going to prove now that if these properties take place 
at bare level they remain true all along the flow. This implies, in particular, that 

all rf; ' n ' are zero, due to the non-existence of any incoming time in that case. 

As explained before, the NPRG equations for any vertex /^ n ' m ' can be rep- 
resented diagramatically by 1-loop diagrams, where vertices and propagators are 
read from Ft + ASk- To those diagrams one must apply the operator dk in order to 
obtain dkrj,"' m \ We represent terms of these equations before the application of 
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the dk operator. Each one of these terms will consist in a number of F^'^ vertices 
joined together by an internal loop of propagators. 

Assuming that the property we want to prove is true at a RG scale fco, we 
will check that it continues to be valid for a scale ko — 8k, which implies the 
result given that by hypothesis it is valid at the bare level. Consider then a typical 

diagram contributing to 9fc.T^ n ' m ' (ti, . . . , t n , ti, ■ ■ ■ ,i m ) as shown in Fig. [14] where 
we have chosen to explicit the time dependence in each vertex. Suppose now that 

dkr[_ n ' m \ti, . . . , t n , ii, . . . , t m ) has a non zero contribution with a given tk strictly 
larger than all ty ■ Let us consider the largest such time tk- Its corresponding 
leg is attached to one of the vertices in the loop, and by hypothesis, at k = fco, 
the largest time in this vertex must correspond to an outgoing leg. However, this 
leg cannot be an external leg, because we have assumed that t k is the largest 
external time. Accordingly, this time (wich we will call t! , see Fig. I14p must be in 
an outgoing internal line. This line arrives to another vertex at a time t" > t! > tk 
that consequently is strictly larger than any external time, and that is associated 
to an incoming internal line. Again, by hypothesis, there must be an outgoing 
leg with a time t'" > t" but this leg cannot be an external leg (because t" > tk)- 
Consequently t'" must correspond to the other internal line that must consequently 
be outgoing. Following the internal (one loop) line and repeating the reasoning for 
each vertex, one returns to the initial vertex at a time strictly larger than the 
supposedly largest time. This is absurd, proving the desired property. 

Now, in the particular case of PA this can be further generalized. In fact, in PA 
we will have that for any /£™' m ) the first three times in the vertex must correspond 
to incoming <f> legs (when possible, that is, when n > 2). This follows in a similar 
way: the property is trivial in PA at the bare level, and we will assume by hypothesis 
that it is true at the scale fco- Again, consider a typical diagram contributing to 

9i-T^ n ' m ' ) (ti, . . . ,t n ,ti, . . . ,t m )- By the property shown in appendix [C] we know 
that the two smallest external times ti and ti correspond to two incoming legs in 
the same vertex in the 1PI loop. 

Let us first assume that the third time tz is also attached to this same vertex. 
By hypothesis, at fco, it must correspond to an incoming leg. If it is an external 
leg the result follows. If it is an internal loop leg, then it must be an outgoing leg 
from another vertex in the loop. But again by hypothesis at fco, there must be an 
external incoming leg in this other vertex with a corresponding time smaller than 
t3, which is absurd (unless n < 3, which would correspond to a particular case). 

The other possibility is that the third smallest time tz is attached to another 
vertex. If this time corresponds to an outgoing leg, we will have that in order for this 
vertex not to have an incoming external leg (which would by hypothesis correspond 

to a time smaller than ts) then this must be a rf, 2 ' 2 ^ vertex with both its incoming 
lines being internal loop lines. But then the same argument as before applies, and 
we can find an smaller time in yet another adjacent vertex. This completes the 
proof. 



E NPRG equations for the simplest vertices 



The first consequence to be taken from the causality results in the previous ap- 
pendix is the non-renormalization of the 2-point function i^' . As can be easily 

seen, the tadpole diagram (that in this case includes rj; 2,2 \ see Fig. I15|) is zero due 
to these causal properties: in order to exist, one of its outgoing legs must be associ- 
ated with a time to smaller than one of the incoming legs, which is impossible in PA. 
The remaining possible diagram would include a vertex and is accordingly 

zero. Notice also that the other 2-point functions rj, ' 2 ^ and r^ 2 ' ^ remain equal to 
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dk 




Fig. 15 Tadpole diagram which would contribute to d k r^'^ forbidden by Ito. 



zero in the shifted version (Eq. [9} of PA. In the first case because of property (1151) . 
and in the second case due to causality (see apendix |D|) . 

(2 1) f2 2) 

We now explicitly solve the flow equations for the functions rt ' and r£ . 

(2 2) 

The flow equation for ' reads 

d k r^ a) (pi,p 2 ,P3,P4) = / d k Rk{q)G 2 k {q)G k (pi+P2~q) 



x r { h 2a) (p 1 ,p 2 ,-q, -pi -P2 + g)r fc (2 ' 2) (g,pi +P2- q,-P3,-pA) (80) 

which is a closed form equation, given that the other possible term is a tadpole 
with a _r( 3,3 ) vertex evaluated at an uncausal configuration. This equation is not 

only closed but, moreover, we can show that it implies that r < j^ 2 ' 2 \p\,p2,Pd,,PA) is 
a function of p\ + p2 only. This is trivially the case in the microscopic initial scale, 

where r^' 2 \pi,p2,P3,p4) = 4A. Now, if this property is true for all k > ko, we can 
define 

lk {pi+P2) = -r^l' 2) {p 1 ,p2,p 3 ,P4,) (81) 

and at this scale the flow equation becomes 

d k r {2 ' 2) \ {pi,p2,p 3 ,Pi) = 16Zfc (pi +p 2 ) / d k R k (q)\ G 2 ko (q)G ko (pi + p 2 - q) 
i k = ko j q ifc=fe 

(82) 

showing by iteration that rji 2 ' (j>i,P2,P3,Pa) remains a function of p\ +P2 for any 
k. We then generalize l k (pi + P2) as in (|81[) to any scale k. This definition allows 
one to arrive to the simpler flow equation 



d k l k (p) = K{p) / d k R k (q)Gi(q)G k (p - q) (83) 

J q 

which can be re-written as 

d k {l k \p)) = 28 k f G k (q)G k (p - q) (84) 

J q 

with solution 

hip) = ? ^ (85) 

1 + 2A J q [G k {q)G k (p -q)- G A (q)G A (p - q)) 

Here the dependence on the scale k of the propagator G k enters only through the 
regulator function R k . This is completely equivalent to the results in the main text 
for k = 0. The general NPRG structure is very reminiscent of the case of the O(N) 
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scalar field theory in the large N limit, where one can also show I51U49| that the 
NPRG equations can be solved exactly, on similar lines. 

As for the three point function, one can start in an analogous way from 

its flow equation 



'kr^ ,1) (pi,p2,p 3 ) = / d k R k (q)Gl(q)G k (pi +p 2 -q) 

J q 



f2 '2,') (2 11 

x Jfe ' (piiPa, -q, -Pi -P2- q)ri ' '(q,pi +p 2 - q, -pi ~P2j- (86) 

(2 2) 

Given the knowledge of the equation is again closed. Moreover, one can now 

(2 2) 

perform the same type of analysis that has been done for rt ' in order to prove 
that r^ 1] { Pl ,P2,P3) is a function of pi +P2 and define 

g k ( P i + P2 ) = \r { k 2 ' 1 \pi,P2, P3 ) (87) 

whose equation reads 

d k g k (p) = 4g k (p)X k (p) f d k R k (q)G 2 k (q)G k (p - q). (88) 

J q 

Then, by using (|S3f) . one deduces that for all k 

dk (iM) = %M _ g* Wa>(p) = p. (89) 

Vifc(p)/ lk(p) i 2 k (p) 

As gA(p) =21 = 21a(p), one concludes that 

g k (p) = 2l k (p) (90) 

with the propagator defined as the regulated bare propagator, which in the case of 
a continuum theory would read 

G k ( q , u) = Wn ( 91 ) 

q 2 + iu + R k (q) 

All these results are strictly equivalent to those given in Section [3] of the main 
text in the limit k — ► 0. 

Within this formalism we can also recover the results given in Section [4] of the 
main text with respect to the threshold of the active-to-absorbing transition, by 
following similar lines as before. 

Notice tough that this NPRG approach is not sufficient to construct a a- 
expansion with closed equations. Indeed, to be able to do that we should have 
closed expressions for any vertex in PA, in order to re-write vertices in BARW-DP 
as an expansion in terms of PA vertices. However, the causal properties which we 

used to ignore the tadpole diagrams in the flow equations for -T^ 1 ' 1 ', -f^ 2 ' 1 ' and 

(2 2) 

' are not general enough to allow us to close the flow equations for vertices 
with a higher number of legs. This is problematic even at first order in a in the 
BARW-PC case, so that we are not able to recover our results for d a within the 
pure NPRG formalism. 
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F Scaling dimension d CT at 1- and 2-loop order 

In this appendix we show how to recover the perturbative results of [8] for d a in 
BARW-PC within our formalism. Let us begin by recalling from the main text the 
equation for a, which reads 

a(p) = j-A^G(q)G(p-q)a(q)l(q) (92) 

Using this expression we can recover the 1-loop and 2-loop results for d a (and 
hence d c ). In order to do so it is convenient to get rid of the bare level dependence 
on A and a by writting some sort of RG flow equation. The easiest way to do this 
is by performing a logarithmic derivative w.r.t. v, which yields 



vd v a(p) = - ±vd v \ a(q)l(q)G(p - q)G(q) 

Jq 

=4iu [ &(q)l(q)G 2 (p-q)G(q) (93) 

J q 

This equation can be compar ed with the corresponding RG equation for l(p), ob- 
tained by differentiating Eq. (JT^J 

vdul(p) = -2i\ P )vd v f G{p - q)G(q) (94) 

J q 

At 1-loop order the (/-dependence in a should be weak. In the IR one expects 
this dependence to be dominated by the external momentum p, given that the 
momentum integral is regular for the values of d we are interested in 

vd„a(p) ~ -A»{p)l(p)vd u f G{p - q)G(q) (95) 

Jq 

which, together with (|94p yields 



Accordingly, and given that 



we have the result 



Hp) (97) 



p( P ) I s 



Hp) ~ ^i 2 (p) (98) 



This behaves, when u, |p| 2 < A 2/(2 ~ d) , as 

Hp) ~ |p| 2(2 " d) (99) 

Given the definition of d a (see Section [5j) we find 

du 3c? 4 2 3c (100) 

With this expression we find that d a changes sign at d c = 4/3 as expected. 

Going now to 2-loop order it is convenient to use the logarithmic derivative of 

a 

vdAo^HP) = -4i / l M^-G 2 ( P - q )G(q) (101) 

Jq Hp) 
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We now introduce in the r.h.s. of this equation the 1-loop result <r(p) ~ l 2 {p), to 
obtain 

vd„\og&(p) = -Aij P^G 2 (p- q)G(q) (102) 

We also need the scaling form for l(p), as given by (|64|) in the main text. At 
this point it is enough in order to obtain d a to restrict to p = (given that 

„ , s d-d„ 

a(p,v)~v 2 ) 

d-d a = -8iv V [ ' / ^ j—i- G(q,o;)G 2 (-q, v - u) (103) 

i (2) Jq \ W ) 

which can be evaluated to yield 

d a = 2-3e + 31og Qj e 2 (104) 
the known 2- loops result [8], which corresponds to d c =~ 1.1. 
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